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1. Introduction 


Advanced sensors and computer models for the surveillance, detection, identification, and 
tracking of sound-emitting targets in combat are rapidly emerging (1,2). Most new Army 
acoustic systems use unattended microphone sensors to construct small ground-based, beam¬ 
forming arrays to determine line-of-bearing angles, e.g., the remote netted acoustic detection 
system (3). At the same time, point-to-point transmission of acoustic signals is greatly 
influenced by variations in the wind flow, temperature field, and turbulence over varying 
topography and surface energy budgets (4-6). Thus, the Anny is looking to implement the best 
possible acoustic propagation models, particularly those that will incorporate the impacts of the 
environment (e.g., the atmosphere, turbulence, and terrain) on acoustic signatures. 


A key element in detennining point-to-point acoustic transmission is modeling the variation of 
the effective speed of sound through the lower atmosphere. The effective speed of sound (c e /J), 
which accounts for the bearing of the sound source from a given receiver and the compass 
bearing of the mean wind, is readily derived from measured or modeled estimates of air 
temperature, relative humidity, and wind velocity. It is useful to define the effective speed of 
sound (c e ff) from the following expression given by Noble (7) and Osteshev (5): 


c eff ~ c o + 11 COS(0 W - 71-6 R ) , 


in which c Q = J— 


^ ? i? = 8314.32 J mol 1 iT 1 is the universal gas constant, T is air 
M 


temperature in Kelvin, M is molecular mass, y s = — is the ratio of specific heats, u is the mean 

c v 

of the horizontal wind, 6 W is the bearing of the wind from North, Or is the bearing of the receiver 
from the source, and u cos (0 W - n- 0 R ) is the component of the sound speed along the direction 
of propagation from the source to the receiver. As discussed in Osteshev (5), the effective sound 
speed in equation 1 is valid only for nearly horizontal propagation angles (5). In addition, Wong 


and Embleton (9) have deduced the ratio of specific heats and molar mass 



as a function of 


temperature and humidity in the following form: 


= 0.04833 + (R h - 0.023U , (2) 

M 

in which A, = 9.2 x 10 5 + 5.5 x 10 6 1 + 4.25 x 10 7 r, f?/„ is relative humidity (in the range 0.0 to 
1.0), and t is air temperature (°C). Then, for applications in outdoor acoustics, sound speed 
differences across vertical layers will generally cause acoustic waves to be refracted upward if 
the effective sound speed decreases with height, and they will be refracted downward if sound 
speed increases with height. 
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Monin-Obukhov (M-O) similarity models (10-14) have been very well suited for such acoustic 
applications. For example, the energy budget and meteorological models implemented in the 
Acoustic Battlefield Aid (ABFA) ( 15,16) are relatively simple and quite robust for generating 
low atmosphere profiles of effective sound speed, especially for daytime convective 
atmospheres. However, for the nighttime case, the model in ABFA assumes that the net surface 
radiation is constant with time, reduced only in the presence of clouds, and is balanced entirely by 
the surface heat flux. Therefore, in this report, we are interested in testing the effectiveness (for 
this application) of an earlier quartic solution for the surface temperature (17,18), i.e., an 
alternate energy budget model, which explicitly computes the net (long wave) radiative flux and 
the turbulent fluxes of heat and moisture close to the ground. Then, the primary energy budget 
model results can be used to estimate profiles for wind speed, air temperature, and relative 
humidity, from which, effective sound speed information can be generated. 1 


2. Model Equations 


2.1 A Quartic Solution for the Surface Temperature 

The alternate surface energy budget model described in this section is a semi-empirical, 
equilibrium model (17,18) that explicitly computes each energy balance component, except the 
flux of energy into and out from the ground, which is computed simply as a percentage, 10% 
daytime and 5% nighttime, of the computed net radiation (21). The model's primary input are 
reference level (1 to 2 m) values for wind speed, air temperature, relative humidity, atmospheric 
pressure, surface thermal emissivity, soil water content, and cloud cover information. (For 
daytime calculations, the day of interest, longitude, and latitude of the site are required to 
determine the solar declination, hour, and zenith angles.) The model's primary output are 
calculations of the net radiation, the effective surface temperature, and the soil, sensible, and 
latent heat fluxes. 

In accordance with Carson (22), it is assumed that the energy budget close to the ground (as 
shown in Figure 1) can be written as 

R net = H + L'E + G , (1) 

in which R net is the net radiative flux, H is the turbulent sensible heat flux, L E is the latent heat 
flux attributable to surface evaporation, and G is the flux of heat into or out of the soil. Positive 
values of R net and G are chosen to represent fluxes directed downward while positive H and L E 
values represent fluxes directed away from the air-soil interface. 


'Note that Swanson and Reichard (19) and Swanson (20) have also reported on the use of M-O similarity-type formulations to 
generate low atmosphere sound speed profiles. However, they based their estimates of the M-O scaling constants on data 
retrieved from a ground-based micrometeorological instrument rather than from energy budget principles, as discussed in this 
report. Nevertheless, their goal was very similar, i.e., to incorporate similarity-based estimates of effective sound speed within 
probability of detection (sensor performance) algorithms for unattended acoustic sensor arrays. 
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WET GROUND J 

Figure 1. The radiation and energy budget close to the ground (illustrated by 
F.V. Hansen, 1993). 

The formulations in the model used to compute the incoming short wave solar radiation (Ryi) for 
daytime cloudless skies are patterned after Meyers and Dale (23) and are augmented as needed 
with empirical results by Haurwitz (24) to account for cloudy skies, i.e., 

Rsi = I = IoT rT pTwT ACOSZ , (2) 

in which I 0 is the extra-terrestrial flux density at the top of the atmosphere on a surface nonnal to 
the incident radiation, Z is the solar zenith angle, and 7) are the transmission coefficients for 
Rayleigh scattering (R), absorption by permanent gases (P), water vapor (W), and absorption and 
scattering of aerosols (A). The empirically derived transmission coefficient for cloudy skies (24) 
is defined by the ratio of insolation with partly or completely covered sky to insolation of 
cloudless skies. 

Gates (25) gives the expression for downward long wave radiation (Ru) for clear skies in the 
model, i.e., 

Rif = - 170.9 +1.195 cr Ty , (3) 

2 —8 —2 —4 

in which R/+ is in watts per meter , <7= 5.6697 x 10 W m K , and T r the reference level (1 to 
2 m) temperature in Kelvin. From Paltridge and Platt (26) an addition to equation 3 is made to 
account for clouds, giving the total expression as 

Rif = -170.9+ 1.194 crT ; 4 +0.3 £T c crr c 4 (cc) , (4) 

in which e c is the long wave flux emissivity of the cloud base, T c is the temperature of the cloud 
base in Kelvin, and cc is the fractional cloud cover. Having an estimate of the cloud height, one 
can approximate T c , assuming an average of the dry and moist adiabatic lapse rate. 
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The upward long wave radiative flux (Rit) is computed by a fonnulation from Yamada (27), i.e., 

R L 't = £CTT^ + (l-£)R L i , (5) 

To evaluate equation 5, one requires an estimate of the effective ground or surface temperature, 
T g . We achieve this through the application of M-0 similarity fonnulations (10-14) and the 
flux-gradient hypothesis as advocated by Dyer (12) and Webb (14). 


The similarity-based expressions for the sensible and latent heat fluxes are defined by the 
MARIAH algorithm (28) and are expressed as 


H = -c „ pu* 9* = — 


c p pk 2 V r (T r -T g ) 

<t>m 0h ln " d ln " d 


z 0 


z h 


and 


( 6 ) 


L'E = —V p u * q* = - 


L' pk 2 

0m 0q in 


V r (q r qg) 


z — d , z-d 
-ln- 


z 0 


z h 


( 7 ) 


in which c p is specific heat; p is air density; k is von Karman’s constant; V,- is the reference height 
wind speed; L ’ is the latent heat of vaporization; q, and q g are the reference height and effective 
surface values for specific humidity, respectively; w* is the friction velocity (m/s); ft is the 
potential temperature scaling constant (K); q lf is the specific humidity scaling constant (kg/kg); 


z is height above ground in meters; d is displacement height in meters (—70% of the principal 
roughness element); z 0 is the aerodynamic (surface) roughness height in meters (-14% of the 
principal roughness element); z/, is the thermal roughness height in meters (-13% of the 
aerodynamic roughness); and (j) m , fa, and fa are the dimensionless wind shear, dimensionless 
temperature lapse rate, and dimensionless humidity lapse rate, respectively. Here, the 
dimensionless tenn for humidity is assumed to be equal to the dimensionless temperature lapse, 
i.e., fa = fa, even though field observations have shown that this is not always correct. Also, for 
stable and very stable conditions at night, Mahrt (29) recommends setting z 0 = z/,. Equations 6 
and 7 can be evaluated in a straightforward manner, given expressions for the dimensionless 
shear and lapse rate terms. 


Following the recommendations of Hogstrom (30), we use 

(j) m =[l-19(z/Z)]“ 1/4 and fa = 0.95[l-11.6(z/Z)]“ 1/2 , (8) 

for daytime unstable atmospheric conditions, and we use 

(j) m = 1 + 5.3 (z/L) and (j)^ = 1 + 8 (z/L) (9) 
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for nighttime stable atmospheric conditions”, where Busch (11) defines the M-0 scaling ratio as 


L 0 V u i 


( 10 ) 


in which 0 V - T,( 1 + 0.61c/, ) is the virtual potential temperature and 0 V * = 0* + 0.61 T,q* is the 
virtual potential temperature scaling constant (q* is the specific humidity scaling constant). 
Here, reference level specific humidity, q r , is detennined as a function of relative humidity (/?/,), 
atmospheric pressure (P r ), and air temperature ( T r ) with the formulation from Rogers (31), i.e., 


3.8 R h 

q r = —— cx p 


5.44xl0 : 


V 


1 


1 


A 


2.73.15 T, 


(ID 


r J J 


In contrast, surface specific humidity, q g , is determined as a function of the following parameters 
(32): reference level specific humidity (q r ), saturation specific humidity at the surface (q s ) (i.e., 
substitute into equation I I R/, = 100% and T r = T g ), the water capacity of the field (Of), and 
finally, soil water content (w s ). As a result, surface specific humidity can be expressed as 


in which 


q g =(l- S) q r +5 q s , 


( 12 ) 


8 = 0.25 


1 - cos 


f „ W 1 
n w r 


0 


f )) 


for w s <6f 


(13) 


and 8= 1 for w s > Of . 

Next, we substitute equations 6 and 7 for H and L ’E into equation 1 and define the soil heat flux 
(G) as G = (jR iu -, (where R net = RlI - Rl t + Rsi ~ R.s'l and Rsl = c/.Rsi)- Then after some 
rearranging, we solve equation 1 for T g , which gives the quartic form for the equilibrium model 
as 


(1 - pif <7 Tg — (l-a)R sl -£R Ll )+ T, 


c„ p k z V r 


, , . z-d . z-d 

0m 0h l n - In- 


z h J 


c„ p k l V r T r 


_ L' pk 1 V r (q r -q g ) 

, , , z-d , z-d , , , z-d , z-d 

0m 0h In- In- 0m 0h In- In- 


z h 


z h 


(14) 


in which a is the surface short wave reflectivity, commonly called “albedo,” and (3 is the 
percentage (e.g., 10% daytime and 5% nighttime) of the computed net radiation used to calculate 
the soil (ground storage) heat flux. Naturally, at night, Rsi = 0. The quartic equation is resolved 
numerically via a Newton-Raphson iterative scheme. 


'Alternately, Webb (14) suggests = (/> h = 1 + 5 (z/L) for nighttime stable conditions. 
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2.2 The Energy Budget Model Implemented in ABFA 

Wilson (75) and Wilson and Szeto (16) report on the development of the prototype model ABFA 
for the analysis of atmospheric refraction and turbulence effects on acoustical detection and 
tracking systems. Among the meteorological input capabilities embedded in ABFA is a surface 
energy budget model for generating low atmosphere profdes of effective sound speed. The ABFA 
energy budget model is based on the empirical fonnulations of Burridge and Gadd (33), as 
discussed in Stull’s (21) textbook (Chapter 7). For the nighttime case, the net long wave radiation is 
given by the following expression: 

^lt-^zl = (0-08pcp)( 1 -0-0.1^-0.3 <t m ~0.6a L ) , (15) 

in which a H , &m, and (Jl are the cloud cover fractions for high, middle, and low clouds, 
respectively. (In contrast, cloud cover is accounted for in the parameterization for R lI shown in 
equation 4.) Clearly, this is a much simpler formulation than the solution of the quartic. 

In ABFA, the net surface radiation is balanced entirely by the surface heat flux, Q , so that the 
M-0 temperature scaling constant is calculated as 


6* 


and the M-0 scaling length is derived as 


Q 

pc pU* 


(16) 


L = - 




@ pc p 


(17) 


where the friction velocity can be determined from the reference height wind speed (Vr) and 
values for roughness length (z 0 ) and displacement (d), i.e., 


II* — 


kVy 


^, H ln 


z- d 


(18) 


\ ZO J 

Thus, the ABFA energy budget model output is used to estimate the similarity-based, log-linear 
profdes for wind speed and temperature extended upward through the surface layer, as discussed 
next. 


2.3 The Similarity-Based Profiles 

Given values for the M-0 scaling constants from energy budget considerations, the formulations 
of Benoit (34) and Nickerson and Smiley (35) are applied to determine the similarity-based 
profdes for wind velocity, temperature, and specific humidity as follows: 


V(z) 


V r + 


U* 

k 


- 

f \ 

In 

z 


\ z r v 


■V n 






Vr, 



(19) 
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and 


n \ 0 * 

z) = T r + — 
k 


q{ z ) = 9r + X 
k 


In 


'z' 


In 


vA- y 


'z' 


Vh 


f z^ 




Vh 


\ z r J 


■V a 


v-w 


^ z^ 


v-w 


Wc, 


viy 


^ z > 

viy 


( 20 ) 


( 21 ) 


where, for the nighttime case, the diabatic influence functions for momentum and heat, y/ m and 
y/h respectively, are given by Webb (14), i.e., 


V, 


m 


V,=-5j • 


( 22 ) 


Here, we assume that y/ q = y/ h . 


For extended profdes of relative humidity, we use equation 11 and the following expression for 
atmospheric pressure given by Campbell and Norman (36), i.e., 

Mz-z r ) A 


P(z)=P, exp 


8200. 


(23) 


which yields 


= 100.0 W> exp| 

J.o 


5.44x 10 : 


V 


1 


1 




273.15 T(z) 


(24) 


JJ 


3. Model Results 


3,1 A Hypothetical Case 

A hypothetical nighttime case is outlined in Table 1. Input parameter values were chosen to 
represent low wind speeds and cool temperatures over a relatively smooth surface. For this 
example, the primary output parameters from the alternate energy budget model (described 
previously) are given in Table 2. First, it is shown that fairly good closure exists for this 
calculation, i.e., R net - G is approximately equal to H + V E. This result helps to validate the 

model for partitioning the net radiation into its surface flux components. The small error (i.e., 

—2 

A = 6.83 Wm ) may be attributed to the approximation for soil heat flux. Second, the estimate 
for surface temperature, i.e., T g = T, - 2.25 K, appears to be quite reasonable for moderately 
stable conditions (i.e., L = 35 m). Third, the M-0 scaling constants appear to be valid and useful, 
so that applying equations 19 through 24 provides the nighttime profiles for wind velocity, air 
temperature, and relative humidity for the layer z < 60 m (see Figure 2). 
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Table 1. Energy budget model parameters 


Parameter 

Symbol 

Value 

Reference atmospheric pressure 

Pr 

1000 mbar 

Reference relative humidity 

R h 

40% 

Reference air temperature 

T r 

283.15 K 

Reference wind velocity 

v r 

3.0 ms -1 

Reference height above ground 

z r 

2.0 m 

Principal roughness element or canopy 
height 

h 

0.50 m 

Surface (long wave) emissivity 

8 

0.98 

Soil water content 

W 
vv s 

0.13 m 3 nf 3 


Table 2. Energy budget model results 


Parameter 

Symbol 

Value 

Net radiative flux 

^net 

-86.62 Wm - - 

Soil (ground storage) heat flux 

G 

-4.33 Wm -2 

Sensible heat flux 

H 

-114.50 Wm -2 

Evaporative (latent) heat flux 

L'E 

+ 25.38 Wm -2 

Friction velocity 

u * 

0.35 ms -1 

Temperature scaling constant 

9* 

0.27 °K 

Flumidity scaling constant 

q* 

-2.35E-05 kgkg -1 

Monin-Obukhov scaling length 

L 

34.64 m 

Surface air temperature 

T g 

280.90 °K 


60- 

a) 

60- 

b) 

60- 

c) 

50- 


50- 


50- 


40- 


40- 


40- 


Height (m) 

a 

_,_i_ 

/ 

30- 


30- 

\ 

20- 

/ 

20- 

/ 

20- 

\ 

10- 

J 

10- 

/ 

10- 

\ 

0- 


0- 


0- 



2 4 6 10 12 14 30 35 40 

VMnd velocity (ms -1 ) Air temperature (C) Relative humidity (%) 

Figure 2. Typical similarity-based profiles for the nighttime case derived from the 
energy budget model: a) wind velocity (in units ms -1 ), b) air temperature 
(in units C), and c) relative humidity (in percent). 
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Next, we apply equations 1 and 2 and the data shown before to derive the extended profiles for the 
effective speed of sound (see Figure 3). In Figure 3, the effective sound speed is determined for 
upwind propagation (left), i.e., the bearing of the receiver from the source and the direction of 
the wind from north are the same, and for downwind propagation (right), i.e., the bearing angles 
are 180° opposite. An asymmetry exists in the graph so that the profile for downwind 
propagation (on the right-hand side) is increasing with height above ground level, while the 
curve for upwind propagation (on the left side) is decreasing (slightly) with height. The effect 
that such profile variations in sound speed will have on short range acoustic attenuation is 
discussed next. 



Figure 3. Low atmosphere profiles of the effective speed of 
sound derived from the energy budget model for 
upwind propagation (dashed) and downwind 
propagation (solid). 

3.2 Approximation of Short Range Acoustic Attenuation 

To briefly examine short range acoustic attenuation at night, we use the low atmosphere profiles 
of wind speed, temperature, and relative humidity (shown before) as input to a flat earth, non- 
turbulent acoustic propagation model called the Windows (version) Scanning Fast Field Program 
(WSCAFFIP). WSCAFFIP is a numerical code developed for assessing enviromnental effects on 
short range acoustic attenuation (7,38). WSCAFFIP detennines acoustic attenuation as relative 
sound pressure loss with range and azimuth for a given frequency and source-to-receiver geometry. 
WSCAFFIP contains propagation algorithms to represent the effects of atmospheric refraction, 
diffraction, absorption, and reflection (ground impedance) on acoustic transmission. Table 3 lists 
the model parameters for an initial approximation of short range acoustic attenuation over an open 
grass-covered (h = 0.5 m) field. Figures 4 and 5 show the WSCAFFIP results corresponding to the 
modeled profiles of effective sound speed generated by the alternate (quartic) model. 
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Table 3. WSCAFFIP model parameters 


Parameter 

Symbol 

Value 

Distance of receiver from source 

X 

1000 m 

Range resolution 

Ax 

10 m 

Bearing of receiver from source 

Or 

90° 

Bearing of the wind from North 


90° (upwind) and 270° (downwind) 

Source height above ground 

h s 

2 m 

Receiver height above ground 


1 m 

Frequency of interest 

f 

100- 1000 Hz 

Ground porosity (tall grass pasture) 


0.675 

Flow resisitivity (tall grass pasture) 

a r 

330 kPa s m 2 



Figure 4. Line plots of acoustic attenuation versus range for 
upwind propagation (dashed) and downwind 
propagation (solid), (i.e., WSCAFFIP model results 
for 300 Hz and 600 Hz, only). 
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As stated earlier, acoustic waves will tend to be refracted downward as the effective sound speed 
increases with height and refracted upward as the effective sound speed decreases with height. 
Therefore, one might expect the amounts of attenuation (in part) to depend on the strength of the 
sound speed profde gradient between the ground and the model top. Because generally, the 
effective sound speed profde for upwind propagation has values decreasing with height, we would 
expect to see greater attenuation with range at each frequency for upwind propagation as a result. 
The line plots of relative sound loss with range in Figure 4 (for 300 Hz and 600 Hz only) and the 
contour graphs of attenuation as a function of range and frequency shown in Figure 5 (for 100 to 
1000 Hz inclusive) show this result satisfactorily. 

In contrast, the behavior of the downwind propagation attenuation curves can be understood in 
terms of “ducting” of acoustic modes. For this example, the sound speed profile forms an 
acoustic duct between the surface and 60 m. At lower frequencies (100 to 200 Hz), acoustic 
waves may be too long for any significant ducting to occur. At 300 Hz, there appear to be one or 
two trapped modes, while at 400 Hz (and greater) several modes appear to be trapped, creating 
the interference pattern. In this way, ducting of acoustic energy across the lower boundary layer 
at night is greatly affected by local meteorological profile structure. 

3.3 Comparison to Observations 

At the start of this project, it was thought feasible to present modeled profiles of effective sound 
speed for the nighttime case in comparison to values of effective speed of sound derived from 
selected micro-meteorological data. As an example, the tower data used to generate the three 
time-height series shown in Figure 6 were collected as part of the 1999 Cooperative Atmosphere- 
Surface Exchange Study (CASES-99) (39). CASES-99 was an extensive field experiment that 
focused on the nocturnal boundary layer (NBL) under mainly clear sky, light near surface wind 
conditions (this included occurrences of weak and intermittent turbulence). The experimental 
period was from October 1 to 31, 1999. The test area for CASES-99 was near Leon, KS 
(50 km) east of Wichita, KS. The test site was chosen (in part) because of its relatively flat 
terrain (average slopes were 0.5 degree). A large number of instruments were deployed in an 
area 4.8 x 3.2 km. These instruments included a heavily instrumented 55-m tower, numerous 
10-m towers, multiple radar, lidar, scintillometers, tethersondes, rawinsondes, and research 
aircraft (see http://www.joss.ucar.edu/cases99). 

The subset of CASES-99 data shown in Figure 6 is from the first intensive observing period 
(IOP) on 5 October 1999 and includes the measurements of wind velocity, air temperature, 
relative humidity (5 min averaged) taken on the 55-m main tower approximately 5 km southeast 
of Leon, KS (latitude: 37 deg 38.88’ N, longitude: 96 deg 44.14’ W, elevation 433.7 m above sea 
level) in an open field of mostly dry grasses, approximately 0.15 m in height (on average). 
Surface roughness height for the site was determined to be about 2 cm and thennal emissivity 
was estimated at 0.98 (40). 
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However, we experienced some difficulties in the application of traditional M-0 profile theory 
for calculating nighttime turbulence structure from temperature, humidity, and wind speed 
difference data (see equations 6, 7, and 14). As an example, Figure 7 shows a time series of the 
modeled surface temperature (77), the 5-m measured temperature (T,), and the 5-m measured 
wind speed ( U r ) for the first nighttime IOP. For this example, reference height is z r = 5 m and 
the canopy height is h = 0.15 m. Certainly, the surface temperature estimates are unrealistically 
low. This “runaway cooling,” according to Mahrt (29), is a known defect in computer models that 
attempt to estimate meteorological gradients near the ground at night (during very stable conditions) 
with M-0 similarity functions. To illustrate the problem further, Figure 8 shows a plot (at t = 

0.125 hr past midnight) of the estimated surface temperature (77), the estimated M-0 scaling length 
(L), and the estimated dimensionless wind shear (<p m ) as a function of the 5-m reference height wind 
speed (U r ), as the wind speed is varied from U r = 1.0 ms -1 to U r = 10.0 ms In addition, the data in 
Table 4 and Figure 9 show a more detail view of the results plotted in Figure 8. The data clearly 
show that the model calculation collapses at U r = 3.6 -3.8 ms -1 , where 77 < 5.0 C,L< 8.0 m, and 
(pm >1.6. Then, for still lower wind speeds, L becomes very small—small enough to cause a 
computational sensitivity, as evidenced in the plot for (p m . Given that both <p n and (ph are needed to 
define two principal components of the quartic, it can be reasoned that when L is small and the (p's 
are large, the model defaults to a kind of radiative equilibrium, where little heat energy (sensible, 
latent, or released from the soil) is available to offset the radiative term. When this occurs, i.e., 
when turbulence collapses, the model results are invalid. 

As a final point, we consider the effect on the model results brought about by varied surface 
roughness (see Figure 10). Recall that for this example, the canopy height is h=0.15 m (short 
grass) and the roughness length is zo = 0.02 m. In Figure 10, when U r =3.6 ms 1 and the canopy 
height < 0.25 m, the model again breaks down. Thus, for whatever the order of wind speed and 
roughness length, there appears be a derivable limit beyond which the calculation should be 
discounted. The arrow on Figure 9 suggests such a limit. 

It is interesting to note that others had previously reported difficulties in the use of similarity- 
based models for the nighttime case (in the calculation of optical turbulence intensity \C^]) 
(41-43). As an example, Figure 11 shows results from the CN2 over-land optical turbulence 
model (42,43), derived from two levels (above ground) of conventional micro-meteorological 
infonnation as input (i.e., wind speed, temperature, and humidity retrieved on fannland, over flat 
barren soil, under mostly clear skies). Although the CN2 model estimates are in fairly good 
agreement with the observations during daytime hours, there are several instances at night when 
the model results show either a consistent bias (about a factor of 2 and mainly underestimated) in 
comparison to the measured turbulence data or the model results fluctuate toward extreme 
(minimum) values. Here, the largest errors occur when there is a computational sensitivity, as 
described before. Average surface layer wind speeds for the nighttime periods shown in 
Figure 11 were approximately 3 to 4 ms -1 . 
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Figure 7. Time series of a) modeled surface temperature ( 77,) and 
the 5-m temperature (T r ) and b) the 5-m measured 
wind speed (U,) for the first nighttime IOP. (The dotted 
line in Figure 7b helps to show the large amount of data 
in this example where U < 2.0 ms -1 .) 
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Figure 8. A plot at t = 0.125 hr of a) the estimated surface temperature 
( T g ), b) the estimated M-O scaling length (Z), and c) the 
estimated dimensionless wind shear (</>„,) as a function of the 
5-m reference height wind speed ( U r ). 


Table 4. Breakdown of M-O similarity models for the nighttime case 


U r (ms ') 

T g (C) 

L (m) 

0m 

3.2 

-8.698 

0.005 

837.346 

3.3 

-8.698 

0.008 

547.225 

3.4 

-8.697 

0.016 

275.226 

3.5 

-8.682 

0.066 

71.424 

3.6 

-8.447 

0.275 

18.173 

3.7 

2.579 

4.914 

1.964 

3.8 

5.140 

8.477 

1.559 

3.9 

6.062 

10.945 

1.433 

4.0 

6.647 

13.245 

1.358 

4.1 

7.072 

15.504 

1.306 

4.2 

7.403 

17.780 

1.267 

4.3 

7.671 

20.096 

1.236 

4.4 

7.896 

22.471 

1.211 

4.5 

8.088 

24.908 

1.190 

4.6 

8.255 

27.426 

1.173 

4.7 

8.401 

30.016 

1.158 

4.8 

8.532 

32.693 

1.145 
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Figure 10. A plot for U r = 3.6 ms -1 of a) the estimated surface 

temperature (71,), b) the estimated M-O scaling length 
(Z), and c) the estimated dimensionless wind shear ($„) 
as a function of the canopy height (in meters). 
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Figure 11. Time series of CN2 model output (dashed circle) compared to 

scintillometer data (solid square), determined from input data at 1 
and 4 m above ground collected over a 450-m path on a) 8-9 July 
1992 and b) 13-14 July 1992 (42,43). 
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4. Summary and Conclusions 


The Army is looking to implement acoustic propagation models that will incorporate the impacts 
of the atmosphere, turbulence, and terrain on acoustic signatures. A key element in estimating 
outdoor acoustic attenuation with range (and frequency) is modeling the variation of the effective 
speed of sound through the lower atmosphere. M-0 similarity models have been very well 
suited for this application, especially for daytime convective atmospheres, e.g., the energy 
budget and meteorological models in ABFA. However, for the nighttime case, we were 
interested in testing the effectiveness of an earlier quartic solution for the surface temperature, 
i.e., an alternate energy budget model, which explicitly computes the net surface radiation and 
the fluxes of heat and moisture close to the ground. 

For a hypothetical nighttime case representing moderately stable atmospheric conditions, the 
alternate (quartic) model perfonned quite well. Model estimates for effective surface 
temperature, net radiation, sensible heat flux, soil heat flux, and latent heat flux were comparable 
to earlier published results ( 17 , 18 ). Also, for this case, energy budget closure was met 
satisfactorily. Thus, the primary energy budget model results could be used to predict profdes for 
wind speed, air temperature, and relative humidity, from which, effective sound speed (and 
ultimately acoustic attenuation) information could be generated. 

However, we found that reasonable values for effective surface temperature were not obtainable 
from the alternate model to simulate the first nighttime IOP of the CASES-99 study. 
Computational difficulties were experienced because of low roughness height and low wind 
speed. We found that as the M-0 scaling length becomes small, the energy budget model breaks 
down and the computation is not viable. Thus, until better solutions can be realized through 
future modeling or experimental works, we recommend that the alternate model be applied only 
when input canopy height and reference level wind speeds can affect values for the M-0 scaling 
length larger than about 13 m. 
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